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Summary. Among the effects predicted by the General Theory of Relativity for 
the orbital motion of a test particle, the post-Newtonian gravitomagnetic Lense- 
Thirring effect is very interesting and, up to now, there is not yet an undisputablc 
direct experimental test of it. To date, the data analysis of the orbits of the exist- 
ing geodetic LAGEOS and LAGEOS II satellites has yielded a test of the Lense- 
Thirring effect with a claimed accuracy of 20%-30%. According to some scientists 
such estimates could be optimistic. Here we wish to discuss the improvements ob- 
tainable in this measurement, in terms of reliability of the evaluation of the system- 
atic error and reduction of its magnitude, due to the new CHAMP and GRACE 
Earth gravity models. 

1 Introduction 

The linearized weak-field and slow-motion approximation of the General 
Theory of Relativity (GTR) is characterized by the condition g^ v ~ ?7 M „ + 
h^y where g^ v is the curved spacetime metric tensor, r^,„ is the Minkowski 
metric tensor of the flat spacetime of Special Relativity and the /i M „ are small 
corrections such that /i^l <C 1. Until now, many of its predictions, for the 
motion of light rays and test masses have been tested, in the Solar System, 
with a variety of techniques to an accuracy level of the order of 0.1% 0- It 
is not so for the gravitomagnetic 1 Lense-Thirring effect due to its extreme 
smallness. It can be thought of as a consequence of a gravitational spin-spin 
coupling. 

1 In the weak field and slow motion approximation of GTR the equations of motion 
of a test particle freely falling in the gravitational field of a central spinning body 
are formally analogous to those governing the motion of an electrically charged 
particle in an electromagnetic field under the action of the velocity-dependent 
Lorentz force. In the gravitational case the role of the magnetic field is played by 
the so called gravitomagnetic field which is generated by the off-diagonal terms 
of the metric goi and whose source is the proper angular momentum J of the 
central body. 
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If we consider the motion of a spinning particle in the gravitational field 
of a central body of mass M and proper angular momentum J, it turns out 
that the spin s of the orbiting particle undergoes a tiny precessional motion 

The most famous experiment devoted to the measurement, among other 
things, of such gravitomagnetic effect in the gravitational field of Earth is the 
Stanford University GP-B mission 0j which should fly at the end of 2003. 

If the whole orbit of a test particle in its geodesic motion around M is 
considered as a sort of giant gyroscope, its orbital angular momentum £ un- 
dergoes the Lense-Thirring precession, so that the longitude of the ascending 
node fi and the argument of pericentre u> of the orbit of the test particle are 
affected by tiny secular precessions f2ur, ^lt |HQ 

• 2GJ . 6GJcoai 

"LT = J, ^LT = r, (1) 

c 2 a 3 (l — e 2 )s c 2 a 3 (l — e 2 )? 

where a, e and i are the semimajor axis, the eccentricity and the inclination, 
respectively, of the orbit, c is the speed of light and G is the Newtonian 
gravitational constant. 

Up to now, the only attempts to detect the Lense-Thirring effect on the 
orbit of test particles in the gravitational field of Earth are due to Ciufolini 
and coworkers who analysed the laser data of the existing LAGEOS and 
LAGEOS II satellites over time spans of some years. The observable is a 
suitable combination of the orbital residuals of the nodes of LAGEOS and 
LAGEOS II and the perigee 2 of LAGEOS II according to an idea exposed in 

m 

6f2 L + Cl 6f2 L 11 + c 2 Su L 11 - 60.2^ LT , Ci ~ 0.295, c 2 - -0.35, (2) 

where the superscripts I and II refer to LAGEOS and LAGEOS II, respec- 
tively. The quantity /xlt is the solved-for least square parameter which is 
in Newtonian mechanics and 1 in GTR. The Lense-Thirring signature, en- 
tirely adsorbed in the residuals of fl and Co because the gravitomagnetic force 
has been purposely set equal to zero in the force models, is a linear trend 
with a slope of 60.2 milliarcseconds per year (mas yr _1 in the following). The 
standard, statistical error is evaluated as 2%. The claimed total accuracy, 
including various sources of systematic errors, is of the order of 20% — 30%. 
The main sources of systematic errors in this experiment are 

— the unavoidable aliasing effect due to the mismodelling in the classical 
secular precessions induced on fl and lj by the even zonal coefficients Ji of 
the multipolar expansion of geopotential 

— the non-gravitational perturbations affecting especially the perigee of LA- 
GEOS II |SJ|5]. Their impact on the proposed measurement is difficult to 
be reliably assessed [TU] 

2 The perigee of LAGEOS is not a good observable because clageos — 0.0045, 
while 6lageos ii = 0.014. 
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It turns out that the mismodelled classical precessions due to the first two 
even zonal harmonics of geopotential J 2 and J4 are the most insidious source 
of error for the Lense-Thirring measurement with LAGEOS and LAGEOS 
II. The combination © is insensitive just to J2 and J4. According to the 
full covariance matrix of the EGM96 gravity model the error due to the 
remaining uncancelled even zonal harmonics amounts to almost 13% |12| . 
However, if the correlations among the even zonal harmonic coefficients are 
neglected and the variance matrix is used in a Root-Sum-Square fashion 3 , 
the error due to the even zonal harmonics of geopotential amounts to 46.6% 
|12| . With this estimate and the evaluations of E] for the impact of the 
non-gravitational perturbations the total error in the LAGEOS-LAGEOS II 
Lense-Thirring experiment would be of the order of 50%. If the sum of the 
absolute values of the individual errors is assumed, an upper bound of 83% 
for the systematic error due to the even zonal harmonics of geopotential is 
obtained; then, the total error in the LAGEOS-LAGEOS II Lense-Thirring 
experiment would become of the order of 100%. This evaluations agree with 
those released in [TH] . 

The originally proposed LAGEOS III/LARES mission J^j consists of 
the launch of a LAGEOS-type satellite-the LARES-with the same orbit of 
LAGEOS except for the inclination i of its orbit, which should be supplemen- 
tary to that of LAGEOS, and the eccentricity e, which should be one order of 
magnitude larger in order to perform other tests of post-Newtonian gravity 
|14l [Ht] . The choice of the particular value of the inclination for LARES is 
motivated by the fact that in this way, by using as observable the sum of 
the nodes of LAGEOS and LARES, it should be possible to cancel out to a 
very high level all the contributions of the even zonal harmonics of geopo- 
tential, which depends on cosi, and add up the Lense-Thirring precessions 
which, instead, are independent of i. The use of the nodes would allow to 
reduce greatly the impact of the non-gravitational perturbations to which 
such Keplerian orbital elements are rather insensitive |§] ■ 

In [IQ an alternative observable based on the combination of the residuals 
of the nodes of LAGEOS, LAGEOS II and LARES and the perigee of LA- 
GEOS II and LARES has been proposed. It would allow to cancel out the first 
four even zonal harmonics so that the error due to the remaining even zonal 
harmonics of geopotential would be rather insensitive both to the unavoid- 
able orbital injection errors in the LARES inclination and to the correlations 
among the even zonal harmonic coefficients. It would amount to 0.02%-0.1% 
only ^JE| (EGM96 full covariance and variance RSS calculations). 

3 Such approach is considered more realistic by some authors I10| because nothing 
assures that the correlations among the even zonal harmonics of the covariance 
matrix of the EGM96 model, which has been obtained during a multidecadal 
time span, would be the same during an arbitrary past or future time span of a 
few years as that used in the LAGEOS-LAGEOS II Lense-Thirring experiment. 



4 Lorenzo Iorio 



In regard to the present status of the LARES project, unfortunately, up 
to now, although its very low cost with respect to other much more complex 
and expensive space-based missions, it has not yet been approved by any 
national space agency or scientific institution. 

2 The impact of the CHAMP and GRACE Earth 
gravity models 

From the previous considerations it could be argued that, in order to have a 
rather precise and reliable estimate of the total systematic error in the mea- 
surement of the Lense-Thirring effect with the existing LAGEOS satellites 
it would be better to reduce the impact of geopotential in the error budget 
and/or discard the perigee of LAGEOS II which is very difficult to handle 
and is a relevant source of uncertainty due to its great sensitivity to many 
non-gravitational perturbations. 

The forthcoming more accurate Earth gravity models from CHAMP JS] 
and, especially, GRACE ^5] will yield an opportunity to realize both these 
goals, at least to a certain extent. 

In order to evaluate quantitatively the opportunities offered by the new 
terrestrial gravity models we have preliminarily used the recently released 
EIGEN2 gravity model |2U]. 

With regard to the combination J5J), it turns out that the systematic 
error due to the even zonal harmonics of geopotential, according to the full 
covariance matrix of EIGEN2 up to degree I = 70, amounts to 7%, while if 
the diagonal part only is adopted it becomes 9% (RSS calculation). Of course, 
even if the LAGEOS and LAGEOS II data had been reprocessed with the 
EIGEN2 model, the problems posed by the correct evaluation of the impact 
of the non-gravitational perturbations on the perigee of LAGEOS II would 
still persist. 

A different approach could be followed by taking the drastic decision of 
canceling out only the first even zonal harmonic of geopotential by discarding 
at all the perigee of LAGEOS II. The hope is that the resulting gravitational 
error is reasonably small so to get a net gain in the error budget thanks to 
the fact that the nodes of LAGEOS and LAGEOS II exhibit a very good 
behavior with respect to the non-gravitational perturbations. Indeed, they 
are far less sensitive to their tricky features than the perigee of LAGEOS 
II. Moreover, they can be easily and accurately measured, so that also the 
formal, statistical error should be reduced. A possible observable is 4 

5f2 L + Cl 5f2 L 11 - 48.2/iLT, ci - 0.546. (3) 

According to the full covariance matrix of EIGEN2 up to degree I = 70, 
the systematic error due to the even zonal harmonics from I ~ 4 to I = 70 

4 A similar approach is presented in |19) . although without quantitative details. 
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amounts to 8.5 mas yr _1 yielding a 17.8% percent error, while if the diagonal 
part only is adopted it becomes 5 22% (RSS calculation). EGM96 would not 
allow to adopt J3J| because its full covariance matrix up to degree I = 70 
yields an error of 47.8% while the error according to its diagonal part only 
amounts even to 104% (RSS calculation). Note also that the combination 
(j2J preserves one of the most important features of the combination (J2J): 
indeed, it allows to cancel out the very insidious 18.6-year tidal perturbation 
which is a I = 2, m = constituent with a period of 18.6 years due to 
the Moon's node and nominal amplitudes of the order of 10 3 mas on the 
nodes of LAGEOS and LAGEOS II [2T]. On the other hand, the impact 
of the non-gravitational perturbations on the combination © over a time 
span of, say, 7 years can be quantified in just 0.1 mas yr™ 1 , yielding a 0.3% 
percent error. The results of Table 2 and Table 3 of ^Hj have been used. It 
is also important to notice that, thanks to the fact that the periods of many 
gravitational and non-gravitational time-dependent perturbations acting on 
the nodes of the LAGEOS satellites are rather short, a reanalysis of the 
LAGEOS and LAGEOS II data over just a few years could be performed. 
This is not so for the combination J5J because some of the gravitational 21 
and non-gravitational [H| perturbations affecting the perigee of LAGEOS II 
have periods of many years. Then, with a little time-consuming reanalysis of 
the nodes only of the existing LAGEOS and LAGEOS II satellites with the 
EIGEN2 data it would at once be possible to obtain a more accurate and 
reliable measurement of the Lense-Thirring effect, avoiding the problem of 
the uncertainties related to the use of the perigee of LAGEOS II. 

Very recently the first preliminary Earth gravity models including some 
data from GRACE have been released; among them the GGM01C model 6 , 
which combines the Center for Space Research (CSR) TEG4 model with data 
from GRACE, seems to be very promising. Indeed, the error due to geopo- 
tential in the combination lf2"|. evaluated by using the variance matrix only 
(RSS calculation), amounts to 2.2% (with an upper bound of 3.1% obtained 
from the sum of the absolute values of the individual terms). Instead, the 
combination would be affected at almost 14% level (RSS calculation), 
with an upper bound of almost 18% from the sum of the absolute values 
of the individual errors. However, it should be pointed out that extensive 
calibration tests have to be still performed for the GGM01C model. 



3 Conclusions 

When more robust and complete terrestrial gravity models from CHAMP 
and GRACE will be available in the near future the combination could 

5 The even zonal harmonics are much more mutually uncorrelated in EIGEN2 
than in EGM96. 

6 It can be retrieved on the WEB at http://www.csr.utexas.edu/grace/gravity/ 
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hopefully allow for a measurement of the Lense-Thirring effect with a total 
systematic error, mainly due to geopotential, of some percent over a time span 
of a few years without the uncertainties related to the evaluation of the impact 
of the non-gravitational perturbations acting upon the perigee of LAGEOS 
II. On the other hand, the obtainable accuracy with the combination (J2J, 
whose error due to geopotential is smaller than that of © , is strongly related 
to improvements in the evaluation of the non-gravitational part of the error 
budget and to the use of time spans of many years. 
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